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Hadron masses are computed from quenched twisted mass lattice QCD for a de¬ 
generate doublet of up and down quarks with the twist angle set to 7r/2, since this 
maximally twisted theory is expected to be free of linear discretization errors. Two 
separate definitions of the twist angle are used, and the hadron masses for these two 
cases are compared. The flavor breaking, that can arise due to twisting, is discussed 
in the context of mass splittings within the A(1232) multiplet. 


I. INTRODUCTION 


Twisted mass lattice quantum chromodynamics (tmLQCD) is an extension of the familiar 
Wilson lattice action [l|, Q, O]. TmLQCD is essentially obtained by introducing a chiral- 
flavor rotation in the mass term of the Wilson action for a quark doublet (see Sec. EH), and it 
offers a number of signihcant advantages. TmLQCD does not suffer from the “exceptional 
conhgurations” which arise in quenched and partially quenched Wilson simulations ^ due 
to the existence of unphysical zero modes. In tmLQCD, 0(a) improvement (where a is the 
lattice spacing) can be implemented in a convenient manner. In particular at maximal twist 
0(a) improvement is often automatic that is, without the addition of a clover term j^. 
The chiral nature of the rotation leads to parity mixings, but at maximal twist the mixings 
become mere operator interchanges. In addition, dynamical tmLQCD simulations [3,0] are 
dramatically less expensive than dynamical chiral fermion simulations, and are claimed to 
be comparable to the cost of dynamical s tagg ered fermions P|. 

Nevertheless there are issues |7|, |8|, |l0|, [ill that remain to be understood. Of particular 
concern in this work is the dehnition of maximal twist. The original suggestion of Frezzotti 
and Rossi was to obtain maximal twist by using the critical hopping parameter (kcw) 
from the Wilson action. This procedure (which we refer to as the Wilson Kcw dehnition) was 
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liiQQ- 


However, recent theoretical studies 


used in previously published computations 

ji^ 1^ have suggested that this is generally not precise enough to maintain tmLQCD’s 0(a) 


improvement for small quark masses (for another point of view see Ref. specihc 

alternatives were presented. 

In this work, the effect of using different dehnitions of maximal twist is studied by per¬ 
forming two separate simulations of the quenched hadron spectrum: one using the Wilson 
Kcw dehnition, and the other using a dehnition which arises by requiring that parity mixing, 
induced by the twisted mass term, vanish in a physical basis 

Our spectrum calculation includes pseudoscalar mesons, vector mesons, and baryons with 
spin-parity (1/2)=*" and (3/2)=*= for a degenerate (u,d) quark doublet. Computations were 
carried out at two lattice spacings and four quark masses. 

Flavor breaking induced by the chiral-flavor rotation in tmLQCD should be suppressed 
by powers of the lattice spacing, but need not exactly vanish for a nonzero lattice spacing. 
Such an effect can be observed in principle as a mass splitting between charged and neutral 
pions, though the computation of a neutral pion mass involves diagrams with disconnected 
quark lines and is therefore expensive. Flavor breaking can also arise in a baryon multiplet, 
such as the isospin-3/2 A(1232), which contains states that are not related by u and d 
quark interchange. In this work, we calculate this mass splitting which does not involve any 
disconnected quark line contribution. 

Section mi contains a brief overview of tmLQCD and the determination of the maximal 
twist angle is discussed. Details of the simulations and the results are given in Sec. IIIII The 
discussion of flavor breaking is given separately is Sec. 113 and Sec. 13 contains a summary. 


II. THE ACTION AT MAXIMAL TWIST 


The Lagrangian of continuum QCD for two massless quark flavors is invariant under the 
transformation 

'ip(x) —>• (1) 

where is the third Pauli matrix in flavor space and u is called the twist angle. When 
a degenerate quark mass nig is added, the fermionic part of the Lagrangian is no longer 
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invariant. After the transformation in Eq. (CD it becomes (in Euclidean notation) 

Cf{x) ='rp{x)[]p+m + i^'^^T^'^'ip{x), ( 2 ) 

where is the covariant derivative and the parameters m and p will be called, respectively, 
the standard mass parameter and the twisted mass parameter. They are given by 

m = rrigcosuj, (3) 

fi = rriq sin uj. (4) 

Two choices for the twist angle uj are of particular interest: a; = 0 causes /i to vanish and 
gives the QCD Lagrangian in its usual form, while u = n/2 causes m to vanish and is 
referred to as maximal twist. Note that for the continuum action the twist can not effect 
physics since it can be removed by a held redehnition. 

One immediately sees that vector and axial currents get mixed by a general rotation. For 
example the charged currents. 


1/ (x) = 'u(x)7^d(a;). 

(5) 

A~{x) = u{x)'y^'y 5 d{x), 

(6) 


transform as follows: 



— (x) cos a; — (x) sin a;. 

(7) 

K{x) 

—*• A~(x) cosoj — iV~(x) sinoj. 

(8) 


A tilde has been added to all currents constructed from helds in the twisted basis. Notice that 
the vector and axial currents don’t mix at maximal twist, but their roles get interchanged. 
TmLQCD is obtained when the transformation of Eq. (HD is applied to the mass term of 


the Wilson lattice QCD action 
twisted basis. 


The fermion action used in the present work is in the 


Spii’, D] = a" ^ ^(x) (^M + ^ ^ (y.V. + y.V: - aV:V.)^ ^{x), 


(9) 


where 


Vu'ipix) = Uiy{x)'ip{x + au) —'ip{x), 
W*j^il){x) = il){x) — Ul{x — ai>)ijj{x — ai>). 


( 10 ) 

( 11 ) 
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In the lattice theory, both the mass term and the chiral symmetry breaking Wilson term 
lack invariance under Eq. (HD so, unlike the continuum case of Eq. (ED, the chiral twist can 
not be removed from the action by a held redehnition. Essentially tmLQCD generalizes the 
Wilson fermion action by introducing a relative chiral phase between the mass term and the 
Wilson term. 

The mass parameter M in Eq. m undergoes both additive and multiplicative renormal¬ 
ization, but the twisted mass parameter /i is only renormalized multiplicatively. It is usual 
y, 1^ to divide M into a mass shift and a quark mass excess m which plays the same 
role as the standard mass parameter in Eq. ED- In particular it is the quark mass excess 


which, along with /r, determines the twist angle |2, El- Our quark propagator calculation 
will be done using the fermion action of Eq. (ED in the twisted basis so we dispense with 
this separation into mass shift and quark mass excess and simply quote M as the standard 
mass parameter. 

Various dehnitions of maximal twist are possible, and they will generally differ by lat¬ 
tice spacing effects. The Wilson KcW dehnition of maximal twist is obtained from the 
standard (unimproved) Wilson action (/i = 0) by computing the pseudoscalar meson mass 
as a function of M (equivalently the hopping parameter k), then extrapolating to vanishing 
meson mass. This critical standard mass parameter is used for the mass shift Mgh at 
all /i and simulations at M = M^., fJ- ^ 0 correspond to maximal twist in this dehnition. In 
terms of the hopping parameter, the critical standard mass parameter is 

1 


aMr = 


2kcw 


4. 


( 12 ) 


WUson ... deaniUon ot .n™.l twi.l ha. been need peeviouely in .i.nnlat.on.|iaill 
l14l |. but studies of 0(a) improvement in twisted mass chiral perturbation theory recommend 
against it jhll Q|. Simulations with this dehnition are discussed in the present work as 
well, and are compared to a dehnition of maximal twist arising from considerations of parity 
conservation. 

The tmLQCD action expressed in terms of twisted helds has a parity violating mass term. 
This term may be removed by a held redehnition and then the parity violation is associated 
with the Wilson term. The resulting form of the action is said to be in the physical basis 
j^. The parity conservation dehnition of the twist angle comes from demanding the absence 
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of a mixed pseudoscalar-vector matrix element in the physical basis 





i:(v'-w«)^’no))=o, 


with 

P^[x) = d{x)'y^u{x). 


(13) 


(14) 


Transforming this according to Eq. m and noting that P~^{x) is invariant under the trans¬ 
formation of Eq. (CD , we arrive at 


iEi(h-(2.‘)^’+(0)) 

tana; =-- P 

Ei(4l-(x,l)F+(0)) 


(16) 


where a tilde again denotes currents constructed from helds in the twisted basis. This 
dehnition is essentially independent of Euclidean time t, as will be seen explicitly below. In 
this work, we compute the hadron spectrum for four values of p, with M tuned in each case 
such that oj determined by Eq. is 7r/2. 


III. SIMULATIONS 


A. Maximal twist from Wilson KcW 


In this section, the quenched hadron spectrum is obtained by using the Wilson Kcw 
dehnition of maximal twist. A standard heatbath algorithm with the Wilson plaquette 
action was used to create two quenched gauge held ensembles of 1000 conhgurations each, 
one at /3 = 5.85 and the other at 0 = 6.0, corresponding to lattice spacings near a = 0.123 


m. 


fm and a = 0.1 fm respectively 

Twisted mass quark propagators are computed from the action in Eq. dHD with a periodic 
boundary condition in all directions using a 1-norm quasi-minimal residual algorithm. The 
standard mass parameter for each f3 is set to the value obtained using Kcw from Ref. jl^ . 
and four diherent values for the twisted mass parameter are chosen which span the range of 
quark masses between rris and mg/Q. {rris denotes the strange quark mass.) From the point 
of view of numerical convergence of the propagator calculation, it would be possible to push 
the calculation to smaller quark masses. However, that would probably require even higher 
statistics to get meaningful results so for our hrst pass we have not tried to probe for the 
ultimate small quark mass limit. As it is, we achieve a pseudoscalar-meson to vector-meson 
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mass ratio of about 1/3. Numerical values of the simulation parameters are displayed in 
Table m 

To uncover any finite volume corrections, simulations were run at /d = 5.85 on both 
16^ X 40 and 20^ x 40 lattices. Differences in the mass spectrum were found to be negligible, 
so only results from the larger volume will be presented. Since the smaller volume was near 
(2 fm)^, we chose a comparable physical volume for the (3 = 6.0 simulations. 

The ground state masses of the pion, p meson, (1/2)^ baryons and (3/2)^ baryons are 
obtained from the following operators, dehned in the physical quark basis: 


TT^ 


(16) 

Pf 

= 

(17) 

Pf 

= 

(18) 

P 


(19) 

n 

= ^ahc[dlG'-f^U^]dc, 

(20) 



(21) 


2 1 
~^J^^abc\^aG'^id}^Uc T ~^J^(^abc\^aG'^i^i\dQ^ 

(22) 

A? 

2 1 
~^J^^abc\d(iG'^iUl^dc ~\~ ~^J^(^abc\dnG'^id}^Uci 

(23) 

A- 

^abc\_dg^G'^idl^dci 

(24) 


where C is the charge conjugation matrix, i is a spatial Lorentz index, and sums over color 
indices a, b and c are implied. Since the pi operator tends to produce a somewhat noisy 
correlator in tmLQCD simulations, the operator pi is also considered. 

As is well known, the Aj operators contain a spin 1/2 admixture. At zero momentum, 
the correlator for spatial Lorentz indices is 171 


(o |t (Aj(f,t)Ai( 0 , 0 ))| o) = ^ 3 / 2 ( 7 ) + ^ 7 j 7 iG'i/ 2 (t). (25) 


The spins are easily separated. 


Gi/2{t) = ^Tr[Gji(f)7i7j], (26) 

G^/2{t) = -Tr[G'jj(t)7j7j + Cii(t)], (27) 


where repeated indices are summed. We will use this G 3 / 2 (t) to determine the masses of 
spin 3/2 baryons. 
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To get twisted versions of the operators in Eqs. (HiEl, one simply applies the twist from 
Eq. (^. In practice the simnlation is done by compnting the tmLQCD qnark propagators 
in the twisted basis nsing Eq. © and then applying Eq. m to both ends. The correlators 
then involve the operators nnchanged in the physical basis. 

As discnssed in Ref. baryon two-point correlators conple to both parities even withont 
twisting. For periodic bonndary conditions, the correlator is 


G(t) = S [(1 + ST4)e-"-‘ + (1 - 

s=± 


sG“\t), 


( 28 ) 


where T is the temporal extent of the lattice, and s snms over positive and negative parity. 
This expression is directly applicable to onr tmLQCD simnlations at maximal twist as well. 
We can therefore discnss the gronnd state masses of both parities from the baryon operators 

in Eqs. ffTTM . 

Thronghont this work, all hadron masses are obtained from three-exponential hts to 
the hadron correlators using all time slices except the source. Statistical uncertainties are 
obtained from the bootstrap method with replacement, where the number of bootstrap 
ensembles is three times the number of data points in the original ensemble. 

Our pion mass squared is obtained from the conhguration average of 7r+ and 7r“ and is 
plotted versus p in Fig. ^ Excellent agreement is obtained with the (3 = 5.85 tmLQCD 
results from Ref. j^. A least squares ht of our data to the form 


{am-^Y = ^ + B(an) + C{anY 


(29) 


is also displayed. Notice that the htted curve does not lead to a vanishing pion mass for 
/i = 0, an effect also noted in Ref. [l^. There could be various reasons for this but one’s 
hrst guess might be that it is due to lattice spacing errors in the determination of Kcw and 
of the twist angle. We will revisit this issue, and also consider the pion decay constant, in 
Sec. HUB] in the context of an alternate dehnition of maximal twist. 

Fig. El is a plot of vector meson mass, averaged over and p and over spatial Lorentz 
indices, versus pseudoscalar mass squared, again with tmLQCD results from Ref. for 
comparison. It is interesting to see that the operator from Eq. (HHD has consistently smaller 
uncertainties than the traditional vector operator of Eq. m- It is also noted that a linear 
extrapolation to vanishing pion mass squared produces a vector mass that agrees with the 
experimental p meson mass, at both (3 values. 






Figures ini and 0show the baryon masses obtained from our fits along with positive parity 
results from Ref. [l3| for comparison. In these plots, our nucleon is the conhguration average 
of proton and neutron, while our A is the conhguration average of and 

Extrapolations of the positive parity baryons are within about 10% of the experimental 
masses. The connection of negative parity data to experiment is less clear. 

The 0{a) improvement of tmLQCD should result in better scaling between lattices at 
[3 = 5.85 and [3 = 6.0 than would be obtained from the Wilson action. This was studied 


for the mass spectrum in Ref. 


12j | and discussed for the pion form factor in Ref 


Q 


Our 


present data for the p and nucleon masses have only a very small dependence on lattice 
spacing relative to the statistical uncertainties, as displayed in Figure El 

B. Maximal twist from parity conservation 


The Wilson dehnition of maximal twist used in Sec. Iflf Al relies on Fq. m with the 
numerical values of k,cw taken from Ref. In this section, maximal twist will be dehned 
by tuning using Eq. (ca). To show that the two options are numerically quite different. 
Fig. El displays co as obtained from Eq. m when the standard and twisted mass parameters 
keep the values used in Sec. IIII Al The source is at (Euclidean) time 1 in this plot. Each p 
value is found to correspond to a different twist angle a;, and the angles range from 65° to 
more than 90°. 

In lig ht of cautionary notes against relying on the Wilson Kcw dehnition of maximal 
twist ll|, ll5|, we will now tune uj to 7r/2 according to Eq. (ITHll and recompute the mass 
spectrum. Parameters used in the simulations are listed in Table IH Notice in particular 
that these results use 300 gauge held conhgurations as opposed to the 1000 conhgurations 
used for the Kcw dehnition of maximal twist in Sec. nnn 

Fig. 13 is a plot of the standard mass parameter M tuned such that cn = 7r/2 for each 
twisted mass parameter. The dependence of aM on ap is seen to be quite linear with a 
distinctly nonzero slope. 

The pseudoscalar meson mass squared is plotted for both dehnitions of maximal twist in 
Fig.lHl The least squares ht of each to a quadratic is also shown. It is evident that the parity 
conservation dehnition of maximal twist, i.e. Eq. m, produces a smaller pion mass for 
sufficiently small p, thus demonstrating that the Kcw dehnition did not precisely correspond 










9 


to a; = 7r/2 for this observable. In addition, Fig. |H1 indicates that the parity conservation 
dehnition of maximal twist leads to less curvature in the rri^ versus p extrapolation and a 
smaller intercept for at /i = 0. 

The curvature of as a function of quark mass and the non-zero value at zero quark 
mass are features commonly found in lattice simulations. Higher order lattice spacing errors 
may lead to a residual pion mass even after tuning to maximal twist to remove 0{a) effects. 
Finite volume effects can not be discounted, although our simulation at /? = 5.85 in Sec. 
nnn did not give any indication of such a problem. In addition, the simple extrapolation 
function Eq. ( 121 , although consistent with simulation data, may not be appropriate at very 


small quark masses. Quenched 
of at zero quark mass |19l. 


_chiral perturbation theory suggests a non-analytic behavior 

si- 

The chiral hts of Chen et al. j2^ to a large set of pion masses calculated at very small 
quark masses (using the overlap fermion approach) are very instructive. They show how 
quenched chiral logarithms can conspire to £t simulation data with an essentially linear quark 
mass behavior in with the pion mass dipping down to zero only at extremely small quark 
masses. That paper also contains many references to earlier searches for quenched chiral 
logarithms. 

It may be interesting to look for quenched chiral logarithms with tmLQCD. However 


Ref. 


22j | suggests much more data at quark masses smaller than those considered in this 


work would be needed. If one has to tune the standard mass parameter as a function of the 
twisted mass /i, this would add to the overall cost of such a calculation. 

In addition to the pion mass, our simulations provide a determination of the pion decay 


C0„su„t. The -.d.eet- method of Ref. Q is particularly eouveruent. 


/,= 




(30) 


where the normalization is such that the physical value is f-jr ~ 130 MeV. Fig. El compares 
the results for the two dehnitions of maximal twist. Numerical results from Ref. using 
the Wilson Kcw definition are also displayed. Much like the pion mass squared, the decay 
constant has significantly less curvature in the quark mass when the parity conservation 
dehnition of maximal twist is used, and is consistent with a linear function of as expected 
from chiral perturbation theory. In contrast, data from the Wilson Kcw dehnition are clearly 
not linear and do not approach a realistic value. Fig. ^1 shows that the decay constant 
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is independent of lattice spacing for all quark masses within uncertainties for the parity 
dehnition, while the Kcw dehnition does not scale as well. 

Beyond the pseudoscalar meson, there is very little difference in the computed spectrum 
between the two dehnitions of maximal twist. The p meson, spin-1/2 baryons and spin-3/2 
baryons are displayed in Figs. El and El respectively. Clear signals are obtained for the 
p, nucleon and A hadrons, but data for the negative parity states do not show a convincing 
signal. The p meson and nucleon masses can be seen to be independent of lattice spacing in 
Fig. 1131 over the whole range of quark masses that were investigated. 

IV. FLAVOR BREAKING AND THE A(1232) MULTIPLEX 

For the Wilson lattice action (p = 0) with degenerate quark masses, the up and down 
quark propagators are identical so hadron masses within any isospin multiplet are degenerate. 
In tmLQCD, the up and down propagators (let’s name them U{x,y) and D{x,y)) differ by 
having opposite signs for the p term, so conservation of isospin is not guaranteed. However, 
tmLQCD does satisfy 

U{x,y) = 75 T)^(p,a:) 75 , (31) 

D{x,y) = 75U^(i/,x)75. (32) 

This leads to the equality of any given hadron two-point correlator with its up-down inter¬ 
changed hermitian conjugate. Since masses are obtained from the real parts of correlators, 
we arrive at the following degeneracies: 

and 7r“, 
p'^ and p~, 
proton and neutron, 

A^^ and A“, 

A"*" and A°. 

The difference between charged and neutral pion masses offers a measure of isospin violation, 
though it would require the computation of disconnected quark diagrams. Here, we consider 
a less expensive alternative: mA++.- — m^+,o. 

Figs. ITHl and ITHl show our results for this A splitting in physical units, as obtained respec¬ 
tively from Wilson Kcw rnaximal twist and parity maximal twist. No signihcant difference 
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is found between the two definitions of maximal twist for this observable. Based on the 
better statistics of Fig. (1000 conhgurations rather than only 300), we see evidence that 
mA++.- — m ^+,0 decreases with decreasing lattice spacing, as expected. 


V. SUMMARY 

Twisted mass lattice QCD offers signihcant advantages over the usual Wilson or clover 
actions. First “exceptional conhgurations” are eliminated and secondly there is the possi¬ 
bility, by tuning the mass parameter, of removing 0{a) lattice spacing errors automatically. 
Our results add to the evidence that these advantages are realized in numerical simula¬ 
tions. The spectrum of hadron masses was computed in tmLQCD for two lattice spacings 
and two dehnitions of maximal twist (where the spectrum is expected to be automatically 
0(a) improved). Four values of the quark mass were used, the lightest corresponding to a 
pseudoscalar-meson to vector-meson mass ratio of about 1/3. No problems with convergence 
of the quark propagator calculation were observed so the quark mass could be pushed even 
lower. 

There are different ways to dehne the twist angle and hence determine maximal twist. 
With the parity conservation dehnition of maximal twist, Eq. ina, the pion decay constant 
was found to be independent of lattice spacing for all quark masses, though this was not true 
for the dehnition of maximal twist that relies on the Wilson critical hopping parameter. No 
signihcant lattice spacing dependence in vector meson or nucleon masses was observed over 
the whole range of quark mass for both dehnitions of maximal twist. Additional calculations 
at more lattice spacings are needed before dehnitive statements about scaling can be made. 
However, it is already clear that maximally-twisted tmLQCD is competitive with other 
improved actions and, of course, much better than the unimproved Wilson action. 

Phenomenologically reasonable results were obtained for the pion, p meson, nucleon and 
A(1232). The tensor operator was found to be useful for attaining the p mass with a smaller 
uncertainty. Hints of the negative parity baryons were also obtained. To do better with 
the baryon resonances probably requires a more sophisticated treatment with operators 
optimized for particular states |23l |. 

For a hxed value of the twisted mass parameter /i, the parity conservation dehnition of 
maximal twist produced a smaller pion mass than the Wilson Kcw dehnition. Since the 
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pion mass squared should be minimized at maximal twist, this suggests that the parity 
conservation dehnition is a better dehnition, at least for this observable. In addition, only 
the parity conservation dehnition led to a realistic extrapolation for the pion decay constant. 
Note that implementing this dehnition required tuning of the standard mass parameter for 
every twisted mass value which adds to the overall cost of doing the calculations. For the 
rest of the mass spectrum, our data do not show a statistically signihcant diherence between 
the two dehnitions of maximal twist. 

A lattice artifact of tmLQCD is havor symmetry breaking. This can manifest itself 
as a mass diherence between charged and neutral pions. However, the calculation of the 
neutral pion mass in the presence of havor symmetry breaking is a formidable task since it 
involves disconnected quark line contributions. In this work we considered mass splittings in 
the isospin-3/2 A(1232) multiplet which involved no disconnected quark line contribution. 
From our higher statistics run there is evidence that the havor symmetry breaking ehect is 
present and that it decreases at smaller lattice spacing. 
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TABLE I: Parameters used in the simulations. Parameters in parenthesis denote the one instance 
where two definitions of maximal twist happen to coincide so new data were not generated. 


/3 

^sites 

^configurations 

aM 

afj, 

5.85 

20^ X 40 

1000 

-0.9071 

0.0376 




-0.9071 

0.0188 




-0.9071 

0.01252 




-0.9071 

0.00627 



300 

-0.8965 

0.0376 




(-0.9071) 

(0.0188) 




-0.9110 

0.01252 




-0.9150 

0.00627 

6.0 

20^ X 48 

1000 

-0.8135 

0.030 



-0.8135 

0.015 


-0.8135 

0.010 


-0.8135 

0.005 

300 

-0.8110 

0.030 


-0.8170 

0.015 


-0.8195 

0.010 


-0.8210 0.005 
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TABLE II: The mass spectrum and pion decay constant from our simulations. The vector meson 
masses are obtained from or whichever has the smaller uncertainty. Spin 3/2 

masses are given as averages over pairs within the multiplet. 


p 

aM 


am^ 


amp 

amjsf 


amA+,A0 

5.85 

-0.9071 

0.0376 

0.4341(4) 

0.1144(7) 

0.647(4) 

0.943(3) 

1.079(6) 

1.056(6) 


-0.9071 

0.0188 

0.3060(4) 

0.0995(7) 

0.566(8) 

0.797(5) 

0.976(9) 

0.933(11) 


-0.9071 0.01252 

0.2518(5) 

0.0923(7) 

0.534(11) 

0.741(7) 

0.939(11) 
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FIG. 1: Pseudoscalar meson mass squared as a function of the twisted mass parameter. The 
standard mass parameter was held fixed to the value obtained from the Wilson KcW definition of 
maximal twist, Eq. m- Data labelled by “xLF” are taken from Ref. 
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FIG. 2: Vector meson mass as a function of the pseudoscalar meson mass squared. The standard 
mass parameter was held fixed to the value obtained from the Wilson KcW definition of maximal 
twist, Eq. (O. Data labelled by “yLF” are taken from Ref. Q- 
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P=5.85 P=6.0 



FIG. 3: Spin 1/2 baryon masses (both parities) as functions of the pseudoscalar meson mass 
squared. The standard mass parameter was held fixed to the value obtained from the Wilson KcW 


definition of maximal twist, Eq. m- Data labelled by “xLF” are taken from Ref. 
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FIG. 4: Spin 3/2 baryon masses (both parities) as functions of the pseudoscalar meson mass 
squared. These data are averages of the A"*"*' and A“ masses. The standard mass parameter was 
held fixed to the value obtained from the Wilson KcW definition of maximal twist, Eq. m- Data 


labelled by “xLF” are taken from Ref. 
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(a/r„)" 


FIG. 5: Vector meson and (1/2)"*' baryon masses as functions of lattice spacing squared. Axes are 
in units of tq = 0.5 fm. The calculation was done using the Wilson HcW definition of maximal 
twist, Eq. m- 
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FIG. 6: The twist angle obtained from Eq. m at different lattice time slices while using the 
Wilson KcW definition of maximal twist, Eq. (d. 
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P = 5.85 |3 = 6.0 



FIG. 7: The standard mass parameter as a function of the twisted mass parameter, tuned such 
that the twist angle oj of Eq. (n3|) is 7r/2. 
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FIG. 8: Pseudoscalar meson mass squared as a function of the twisted mass parameter, calculated 
using the parity conservation definition of maximal twist: a; = 7r/2 in Eq. (1151) . Data calculated 
using the Wilson HcW definition. Fig. ^ are shown for comparison. 
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FIG. 9: Pseudoscalar meson decay constant as a function of the pseudoscalar meson mass squared, 
calculated using both the parity conservation definition and the Wilson HcW definition of maximal 


twist. Data labelled by “yLF” are taken from Ref, 
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FIG. 10: Pseudoscalar meson decay constant as a function of lattice spacing squared, calculated 
using both the parity conservation definition and the Wilson KcW definition of maximal twist. Axes 
are in units of tq = 0.5 fm. 
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FIG. 11: Vector meson mass as a function of the pseudoscalar meson mass squared, calculated 
using the parity conservation definition of maximal twist: a; = 7r/2 in Eq. (1151) . Data calculated 
using the Wilson KcW definition, Fig. ^ are shown for comparison. 
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FIG. 12: Spin 1/2 baryon masses (both parities) as functions of the pseudoscalar meson mass 
squared, calculated using the parity conservation definition of maximal twist: cj = 7r/2 in Eq. (|1 .Ij) . 
Data calculated using the Wilson Kcw definition. Fig. ^ are shown for comparison. 
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FIG. 13: Spin 3/2 baryon masses (both parities) as functions of the pseudoscalar meson mass 
squared, calculated using the parity conservation definition of maximal twist: uj = 7r/2 in Eq. (1151) . 
These data are averages of the and A“ masses. Data calculated using the Wilson KcW 

definition, Fig. ^ are shown for comparison. 



FIG. 14: Vector meson and (1/2)''' baryon masses as functions of lattice spacing squared, calculated 
using the parity definition of maximal twist: u = 7r/2 in Eq. (|15l) . Axes are in units of tq = 0.5 
fm. 
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P=5.85 P=6.0 



FIG. 15: Flavor splitting within the A(1232) multiplet as a function of the pseudoscalar meson 
mass squared. The standard mass parameter was held fixed to the value obtained from the Wilson 
Kcw definition of maximal twist, Eq. m- 
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FIG. 16: Flavor splitting within the A(1232) multiplet as a function of the pseudoscalar meson 
mass squared, calculated using the parity conservation definition of maximal twist: cu = 7r/2 in 

Eq. (HHI). 













